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基于 改进 正 余 缀 优化 算法 的 多 国 值 图 像 分 割 
郎 春 博 ， 贾 稚 鸣 :， 邢 致 僧 ， 彭 晓 旭 ， 李 金 夺 ， 康 立 飞 


(东北 林业 大 学 机 电工 程 学 院 , 哈尔滨 150040) 


摘 要 : 针对 多 阅 值 图 像 分 割 方法 计算 量 大 、 分 割 精度 低 的 问题 ,提出 了 基于 改进 正 余 弦 算 法 (improved sine cosine 
algorithm) 的 多 阔 值 图 像 分 割 方法 。 首 先 对 种 群 进行 混沌 初始 化 来 提高 初始 种 群 质 量 ; 其 次 根据 粒子 适应 度 值 的 大 
小 自 适应 地 调整 参数 ; 最 后 引入 反 向 学 习 策 略 并 择优 选取 粒子 。 伯 克利 图 像 和 植物 冠 层 图 像 分 割 实验 的 结果 表明 ， 
该 算法 的 运行 时 间 较 短 ， 而 且 分 割 精度 较 高 ， 具 有 较 强 的 便 棒 性 
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Multi-threshold image segmentation based on improved sine cosine optimization algorithm 


Lang Chunbo, Jia Hemingi, Xing Zhikai, Peng Xiaoxu, Li Jinduo, Kang Lifei 
(College of Mechanical & Electrical Engineering, Northeast Forestry University, Harbin 150040, China) 


Abstract: Aiming at the problem of the computational complexity and low segmentation precision of the multi-threshold 
image segmentation method, this paper proposed an Improved Sine Cosine Algorithm (ISCA) based multi-threshold image 
segmentation method. Firstly, this method used a chaotic initialization technique to improve the quality of initial population. 
Secondly, it introduced an adaptive strategy to adjust the parameters according to the fitness values. Finally, it utilized an 
opposition-based learning strategy, then the better particles were selected. The results of the Berkeley image and the plant 
canopy image segmentation experiments show that this method has a satisfied performance in terms of running time and 
segmentation accuracy. And it has a strong robustness. 
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引言 等 人 0 于 2016 年 提出 的 新 颖 随机 群 优 化 算法 。 该 算法 的 核 
了 心思 想 是 利用 正 余弦 函数 的 震荡 特性 ， 向 外 波动 进行 全 局 探 
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图 像 分 割 是 由 图 像 处 理 到 图 像 分 析 的 关键 步骤 ， 其 目的 索 ， 向 最 优 解 波动 进行 局 部 开发 ， 使 算法 逐步 收敛 。SCA 具 
于 将 图 像 分 割 成 若干 个 独立 的 、 有 具有 特殊 性 质 的 部 分 并 提 ” ”有 易于 实现 ， 收 全 速度 快 ， 求 解 精度 高 等 优点 ， 其 随机 参数 
出 感 兴趣 目标 。 目 前 常用 的 图 像 分 割 方法 有 基于 区 域 的 分 。 与 递减 参数 的 设置 也 使 得 算法 的 探索 与 开发 能 力 得 到 很 好 的 
I 方法 、 基 于 边缘 的 分 割 方法 、 基 于 特定 理论 的 分 割 方法 以 平衡 ， 但 局 部 搜索 能 力 较 弱 的 缺点 使 得 算法 效果 并 不 理想 。 
基于 阔 值 的 分 割 方法 等 上 9。 而 阔 值 分 割 方法 因 其 计算 量 小 、 本 文 提出 一 种 改进 正 余弦 算法 (improved sine cosine 
于 实现 效率 高 、 性 能 稳定 等 优点 被 广泛 应 用 于 很 多 领域 。 algorithm, ISCA ) 并 将 其 应 用 到 多 阔 值 图 像 分 割 领域 。 为 解 
如， 在 交通 图 像 分 析 中 ， 汽 车 牌照 的 识别 与 提取 中， 在 医 。 决 标准 SCA 局 部 搜索 能 力 较 弱 的 问题 , 在 保证 运行 效率 的 同 
学 图 像 分 析 中 ，PET-MR 图 像 判 断 病情 的 恶化 或 好 转 中 ;在 时 提高 分 割 精度 ， 进 行 了 以 下 几 部 分 改进 : a) 混沌 初始 化 种 

学 图 像 分 析 中 ， 研 究 森 林 火 灾 图 像 来 检测 实时 火 情 ， 为 实 。 群 ， 提 高 初始 种 群 质 量 ;b) 采用 一 种 自 适 应 参数 策略 代 蔡 

火灾 的 预测 预报 奠定 坚实 基础 中 。 传 统 的 单 闽 值 分 割 方法 有 的 线性 递减 参数 策略 ， 根 据 粒 子 适 应 度 值 的 大 小 动态 调 
适应 性 并 不 广泛 ， 因 此 一 些 学 者 利用 和 迭代 的 方式 将 其 扩展 至 参数 以 更 好 地 平衡 探索 与 开发 能 力 ;c) 引入 反 向 学 习 策 略 
闵 值 。 加 可 选 粒 子 数量 ， 并 择优 选取 粒子 ， 增 大 随机 性 的 同时 提 
多 阔 值 图 像 分 割 可 以 被 看 成 一 种 优化 问题 。 传 统 多 阔 值 。 种 群 进化 速度 。 实 验 选 取 四 幅 伯 克利 图 像 以 及 两 幅 植 物 冠 
方法 的 本 质 就 是 穷 举 法 ， 其 计算 量 会 随 着 闵 值 个 数 的 增加 时 图 像 ， 对 比 SAPSO、SSA、PSO、ALO 与 HSO 五 种 优化 
指数 增长 ， 实 时 性 较 差 。 为 此 有 学 者 利用 智能 算法 对 分 割 的 法 。 实 验 结果 表明 ，ISCA 增强 了 算法 的 局 部 搜索 能 力 ， 
适应 度 函 数 进行 寻 优 操作 ， 在 一 定 程度 上 解决 了 运行 效率 的 得 求解 精度 得 到 有 效 地 提高 ， 其 运行 时 间 、 分 割 精度 等 指标 
问题 *。 但 算法 本 身 的 局 限 性 ， 寻 优 能 力 并 不 理想 ， 均 优 于 对 比 算法 ， 具 有 较 高 的 鲁 棒 性 ， 为 后 续 的 图 像 提 取 操 
易 早熟 收敛 、 陷 入 局 部 最 优等 状况 仍 需 得 到 有 效 地 改善 。 作 黄 定 了 良好 的 基础 。 

正 余 弦 算 法 (sine cosine algorithm, SCA) 是 由 Mirjalili 
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录用 定稿 郎 春 博 ， 等 : 基于 改进 正 余 蓄 优 化 算法 的 多 阔 值 图 像 分 割 第 37 卷 第 4 期 
me 
1 ”多 阐 值 Otsu 图 像 分 着 正 余弦 函数 间 随 机 切换 。 
Otsu 算法 (最 大 类 间 方 差 法 ) 是 一 种 自 适应 确定 一 值 化 3 ”改进 正 余弦 优化 算法 〈ISCA) 
图 像 分 割 阔 值 的 算法 , 由 日 本 学 者 大 津 于 1979 年 提出 , 该 算 。 ”3.1 混沌 初始 化 
法 易于 实现 且 在 一 定 条 件 下 不 受 亮度 或 对 比 度 影响 ， 因 而 被 初始 种 群 质量 对 SCA 算法 来 说 是 至 关 重 要 的 , 而 一 个 良 
广泛 应 用 于 图 像 处 理 领 域 [3-13。 其 核心 思想 为 : 利用 闵 值 将 。 ”好 的 初始 种 群 可 使 其 在 收敛 速度 、 求 解 精 度 等 方面 得 到 显著 
图 像 分 割 为 背景 、 前 景 两 部 分 ， 使 得 这 两 部 分 类 间 方 差 最 大 ”提升 。 以 往 的 初始 化 操作 是 基于 rand 函数 进行 的 ， 尽管 该 方 
的 阐 值 便 为 最 优 阔 值 。 法 有 着 良好 的 随机 性 ， 但 其 分 布 情况 往往 不 尽 人 意 。 而 混沌 
假设 一 幅 图 像 具 有 MxN 像素 、 工 个 灰 度 级 ， 其 中 灰 度 。 变量 具有 随机 性 、 遍 历 性 、 规 律 性 的 特点 ， 可 以 很 好 地 保持 
级 为 ; 的 点 的 概率 为 群体 多 样 性 [7-181。 这 里 选用 Logistic 这 一 经 典 混沌 映射 进行 
pi = /MN (1) ”初始 化 操作 ， 其 表达 式 如 下 : 
其 中 :ni 表示 灰 度 值 为 i 的 像素 个 数 。 Xn = HN (1 一 为 ) (6) 
设 阔 值 :将 整 幅 图 像 分 为 前 景 和 背景 两 部 分 ， 则 类 间 方 ” 其中:， ws(04] ，x% s(0.0 ，/ 越 大 混沌 性 越 高 ，w=4 时 系统 
差 co 是 : 的 函数 ， 处 于 完全 混沌 状态 。 

ry pe 混沌 初始 化 步骤 如 下 : 

Oe ee a) 利 用 rand 函数 随机 生成 种 群 的 第 一 个 粒子 x(i.:) 。 
Pr St 混沌 初始 化 下 一 个 粒子 , 对 x(i;;) 的 每 个 维度 进行 归 一 
9 Bs 加 

各 气 化 操作 并 赋值 给 xi+1:) ， 使 其 在 [0 内 。 
Le. 3 请 c) 根 据 式 06) 进行 Logistic 映射 ， 生 成 混沌 变量 cx(7)， 
国王 之 名 = 和 mo = Zip 生出 加 2..D， We 
和 分 别 表示 像素 被 分 到 前 景 和 背景 两 类 中 的 概率 ， dj 利用 混沌 变量 ox(j) 将 x(i+17) 映射 回 原始 空间 ， 公 式 
m 和 mi 分别 表 示 两 类 的 平均 灰 度 ，ms 表示 整 幅 图 像 的 平均 。 为 
灰 度 ， 使 ce () 达到 最 大 时 的 7 即 为 最 优 阔 值 。 cs 
将 该 公式 推广 至 多 阔 值 ， 设 阔 值 组 合 为 (ist) ， 此 ee 
时 图 像 被 分 为 上 个 类 别 ， 类 间 方 差 表 示 为 : e) 判 断 是 否 遍 历 所 有 粒子 ， 否 则 转向 b)。 
， 所 3.2 自 适应 参数 策略 
ei 根据 以 上 分 析 不 难看 出 ，SCA 算法 中 参数 起 着 至 关 重 
使 得 oo) 达到 最 大 值 的 阔 值 组 合 (人 7) ”要 的 作用 ， 当 ;>1 时 控制 粒子 向 远离 最 优 解 方向 移动 以 加 强 
即 为 最 佳 阔 值 组 合 。 但 传统 的 多 阔 值 Otsu 方法 其 本 质 就 是 穷 全 局 搜索 能 力 ，i <1 时 控制 粒子 向 着 最 优 解 移动 以 加 强 局 部 
举 法 ， 总 计算 量 近 似 为 0(0)* ， 且 随 着 阔 值 个 数 天 的 增加 时 。 ”开发 能 力 。 但 标准 SCA 算法 中 参数 仅仅 随 着 达 代 次 数 的 增 
指数 级 增长 ， 为 提高 程序 运行 效率 ， 现 采用 ISCA 算法 对 多 ”加 而 线性 减 小 ， 并 没有 考虑 到 每 个 粒子 的 实际 情况 ， 严 重 影 
阀 值 选取 过 程 加 以 优化 。 响 了 种 群 的 进化 速度 以 及 算法 寻 优 能 

玉 此 本 文采 用 了 一 种 自 适 应 参数 调整 策略 ， 该 策略 可 以 

2 正 余弦 优化 算法 《SCA) 针对 每 个 粒子 的 不 同情 况 加 以 赋值 ， 在 保证 全 局 探索 质量 的 

正 余弦 优化 算法 是 由 Mirjalili 等 人 于 2016 年 提出 的 新 的 ”同时 增强 了 SCA 算法 局 部 开发 能 力 09!。 相 较 原来 仅 随 迭 代 
随机 群 优化 算法 ， 其 基本 思想 是 利用 正 余弦 函数 的 震荡 特性 。 ”次数 变化 而 线性 改变 参数 的 策略 而 言 ， 自 适应 参数 策略 普遍 
逐步 收敛 于 最 优 解 ， 向 外 波动 进行 全 局 探索 ， 向 最 优 解 波动 性 更 强 、 适 应 度 更 高 ， 公 式 如 下 : 
进行 局 部 开发 。 本 

利用 SCA 算法 求解 优化 问题 的 思路 为 ， 在 Dp 维 空间 中 a 下 < 加 
初始 化 产生 数量 为 m 的 种 群 , 种 群 中 第 ii=12…m) 个 粒子 的 ey Re Pe 而 
位 置 表示 为 =(zzsxo)， 通 过 计算 粒子 对 应 适应 度 函 数 EN 
值 (3) 并 加 以 比较 ， 得 到 当前 最 优 个 体位 置 rion + Ven) xrand fi> fo 
Pe =(pasPe,… Pw) ， 其 中 粒子 更 新 位 置 的 公式 如 下 : 

pl es 二 其 中 ， /为 第 i 个 粒子 的 适应 度 函 数值 ， 及 为 所 有 粒子 的 平 

”yy +ixcos(s)xlspy -| A205 均 适 应 度 函数 值 ， fi, 为 所 有 适应 度 值 优 于 /的 粒子 其 适应 
其 中 : ne(0,2x) ，ss(02) ，ne(0,1) 为 均匀 随机 数 , 参数 1 线 。” 度 值 的 平均 ， 焉 为 当前 最 优 适 应 度 值 ， ma 的 定义 如 下 
性 递减 ， 其 表达 式 如 下 : Timax + Fimin 

有 fee = (8) 
| - 当 有 < fi 时， 当前 粒子 距离 种 群 中 最 优 粒子 位 置 较 近 ， 
其 中 :< 为 常数 ，; 为 当前 迭代 次 数 ， 心 为 最 大 和 迭代 次 数 。 ”应 赋 给 其 一 个 较 小 的 参数 5 来 加 强 算 法 的 局 部 开发 能 力 ， 训 


由 式 (4) 可 以 看 出 
1 。 正 是 因为 随机 参数 和 线性 
好 地 平衡 了 全 局 探 
下 一 次 移动 的 方向 
入 了 一 定 的 
增 大 或 减 小 当前 移动 方向 所 产生 的 效果 ， 


H SCA 


算法 共有 4 


个 主要 参数 ,分 别 为 i,5,# 


免 错 过 最 优 位 置 ， 当 fh <f<h 时， 当前 粒子 在 种 群 内 质量 


修正 = pe 


递减 参数 的 设置 使 得 算法 很 一般 ， 念 其 
索 与 局 部 开发 能 力 nW91。 参 数 i 决定 了 粒子 ”探索 与 开发 
参数 则 控制 了 移动 的 距离 ， 参 数 * 引 ”位 置 较 不 理想 ， 
策略 ， 赋 给 当前 最 优 解 一 个 随机 的 权 信 ， 以 子 靠近 


参数 ”的 作用 是 在 


参数 主 随 迭代 次 数 增加 而 线性 
能 力 ; 当 f> 太 时 ,当前 粒子 距离 最 


递减 以 较 好 的 平衡 
多 粒子 较 远 ， 


故 赋 给 其 一 


个 较 大 的 参数 以 尽快 向 最 优 粒 
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3.3 反 向 学 习 策 略 
反 向 学 习 策略 是 由 Tizhoosh 于 2005 外 


学 习 方法 ， 


它 的 反 向 解 相 比 有 50% 的 概率 更 远离 最 优 解 P。 


邹 春 博 ， 等 


E 提 出 的 一 种 优化 
产生 的 候选 解 与 
因此 ， 通 过 


该 策略 的 理论 依据 是 : 每 个 随机 


寻找 当 


前 解 的 反 向 解 ， 
秀 的 解 作 为 候选 解 ， 
进化 速度 以 及 算法 的 寻 


其 中 : 避 为 种 群 中 第 i 个 粒子 在 第 j 维度 上 的 反 向 数 ; 


对 二 者 进行 评估 ， 
进而 引导 种 群 癌 最 优 解 各 
能 力 ， 公 式 如 下 : 


=k(aj +b))—xy 


贪 焚 地 选择 较 优 
近 ， 提 高 种 群 


(9) 
a 和 b 


分 别 为 第 j 维度 的 上 界 与 下 界 ; 参数 上 决定 反 向 学 习 的 模型 。 


本 文 依据 一 定 概率 执行 反 向 学 习 操作 ， 虽 
性 的 同时 降 


在 增 大 随机 


氏 算 法 复杂 度 。 


3.4 算法 描述 


混沌 初 
分 布 的 随机 1 
的 适 配 度 ， 


台 化 策略 能 够 生成 高 质量 的 初始 种 群 ， 增 强 粒 子 
生 ;， 自 适应 参数 策略 则 提高 了 粒子 与 其 对 应 参数 
与 仅 随 迁 代 次 数 增 大 而 改变 的 线性 递减 策略 相 


比 , 更 好 地 平衡 ] 
略 增加 了 种 群 多样 性 


SCA 算法 的 探索 与 开发 能 力 ; 反 向 学 习 策 
， 贪 禁地 选取 优秀 个 体 ， 增 强 算法 局 部 


人 开发 能 力 的 同时 加 快 了 种 群 进化 速度 ， 提 高 搜索 效率 ， 其 流 
me 蛙 如 图 1 所 示 。 
OY) 程 如 所 示 
DE 3 混 池 初 始 化 种 妊 ， 计算 各 粒子 
一 2 “给 定 算法 参数 ” | 并。 挝 应 度 舍 f/l0 
©O y l! 
OO 更 新 各 粒子 位 置 。 |< 一 | 自 适应 更 新 参数 /。 | 一 | 。 初始 化 最 优 值 

一 算 各 粒子 2 行 后 下 | 兵 : 
GO 适应 弃 信 f+1) > 执行 反 向 学 习 操 作 
| 否 

ey 更 新 最 优 值 | | 择优 选 出 候选 解 
GN 

人 图 1 改进 正 余 弦 算 法 流程 图 

CC Fig.1 Flow chart of the Improved Sine Cosine Algorithm 
全 改进 正 余弦 算法 (ISCA) 的 基本 步骤 如 下 : 

OO a) 随 机 生成 群体 规模 为 1 初始 种 群 ， 利 用 式 (6) 混沌 

初始 化 每 个 粒子 的 位 置 ， 预 先 设 定 算法 wl Ds 


b) 计 算出 每 个 粒子 适应 度 值 /GD) 


为 有 。 


0o) 根 据 式 《7) 自 适 
dd) 根据 
围 [0.5L-1] 内 。 
6) 将 更 新 后 的 位 置 代入 适应 


f (itl)。 


应 地 调整 每 个 粒子 的 参数 。 
式 (4) 更 新 各 粒子 位 置 ， 同 时 将 其 限制 在 最 大 范 


人 车 rand > p ， 则 根据 式 (9) 进行 反 向 学 习 操作 ， 并 在 当 


前 解 与 反 向 解 中 择优 选 出 候选 解 。 

8 将 候选 解 对 应 适应 度 值 与 最 优 值 比较 ， 若 优 于 天 则 
将 其 更 新 。 

h) 判 断 是 否 满足 迭代 终止 条 件 ， 若 满足 则 退出 循环 ， 否 
则 转向 c)。 


4 ”实验 


实验 背景 
为 了 验证 ISCA 算法 在 多 阔 值 图 像 分 割 中 的 有 效 性 ， 选 


4.1 


结果 及 分 析 


月 未 


取 4 幅 伯克利 经 


图 像 以 及 2 幅 植 物 冠 层 图 像 加 以 分 析 研 究 ， 


应 度 函 数 得 到 新 的 适应 度 值 


基于 改进 正 余弦 优化 算法 的 多 冰 值 图 像 分 割 


并 与 SAPSO、SSA、PSO、ALO、HSO 算法 进行 比较 。 其 中 
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SAPSO 算法 引入 模拟 退火 思想 , 利用 模拟 退火 算法 的 “ 突 跳 ” 


能 力 有 效 地 避免 陷入 局 部 最 优 ， 提 高 了 算法 性 能 ; 
ee 
HSO 和 声 搜索 则 在 近 些 年 被 ) 

程序 运行 时 间 以 及 图 像 分 割 法 性 能 进行 评判 。 


SSA 楼 海 
ALO 蚁 狮 算法 与 
应 用 于 优化 领域 。 本 文通 过 


构 相 似 性 


其 中 对 于 图 像 分 割 精 度 的 分 析 引 入 峰值 信 噪 比 (PSNR) 和 结 


(SSIM) 两 个 指标 ， 其 定义 如 下 ?11; 
PSNR = iomgu| 2 (10) 
MSE = EUGED KEN GD 


其 中 :7(i, 站 和 K(i,j) 分别 表示 大 小 为 MxN 的 原 图 像 和 分 割 
图 像 。 
(2 人 +c)(20, +c,) 
SSIM (x, y) te (oa oom) (12) 
其 中 : .和 分别 为 原 图 像 和 分 割 图 像 平均 灰 度 ; 有 和 吕 
分 别 为 原 图 像 和 分 割 图 像 的 方差 0 为 原 图 像 和 分 割 图 像 的 


协 方差 ; 


cl =6.5025 ， 


c = 58.5225 


植物 冠 层 图 像 采 集 自 东北 林业 大 学 模式 识别 与 智能 系统 


辩 率 为 3000 X4000。 算 法 的 实验 环境 为 


Windows10 系统 ，1.6 GHz 处 理 器 ，8 GB 内 存 ，64 位 操作 系 
统 , 编程 环境 为 MATLAB R2017a。 实验 中 各 参数 设置 如 下 : 


粒子 数目 w=30， 最 大 迭代 次 数 . =500 ， 各 算法 特性 参数 设 
置 如 表 1 所 示 。 
表 1 各 算法 参数 设置 
Table 1 Parameters of each algorithm 
算法 参数 取 值 
nn [0,2] 
ISCA p 0.5 
从 4 
1 -g/lIn(0.1) 
SAPSO 
k 0.95 
SSA a [0,2] 
中 [0.4,0.9] 
We 25:5 
PSO 
商 2 
6 2 
ALO w [2,6] 
PAR 0.3 
HSO HCMR 0.95 
BW 25.5 


4.2 伯克利 经 典 图 像 分 割 实验 


为 了 验证 本 文 优化 算法 在 多 阔 值 图 像 分 割 中 的 普 适 性 ， 


直 图 像 分 割 


克利 图 像 库 中 随机 选取 4 幅 图 片 ， 如 图 2 所 示 。 
差 函 数 作为 算法 的 适应 度 函 数 ， 利 用 ISCA 算法 对 以 上 4 
妈 像 多 阔 值 分 割 的 阔 值 选取 过 程 加 以 优化 。 实 验 中 所 有 图 


以 类 间 


3、4、5。 为 定量 地 比较 各 优化 算法 在 多 
生 能 ， 本 文 对 程序 运行 时 间 以 及 图 像 分 割 


FP 的 


两 方 


田 进 和 


行 了 而 


究 ， 其 中 各 算法 立 值 情况 如 表 2 所 示 。 


E 对 各 


白 
差 
图 
选取 的 阐 值 为 2、 
值 
度 


法 分 割 阅 值 加 以 分 析 : 
几乎 相同 。 当 


当 K=2 时 ， 所 有 算法 的 取 
K=3 时 ，ISCA 算法 搜索 空间 扩大 ， 全 局 搜 


能 力 得 到 提高 ， 而 其 他 算法 性 能 提高 不 明显 。 当 K=4，5 
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时 ，ISCA 算法 的 全 局 搜索 能 力 优 势 得 以 进一步 体现 ， 可 在 
整个 解 空间 内 进行 随机 搜索 操作 ， 有 效 地 避免 了 陷入 局 部 最 
优 。 


其 次 对 各 算法 运行 时 间 加 以 分 析 : 由 表 3 中 数据 可 知 ， 
随 着 阔 值 个 数 的 增加 ， 算 法 复杂 度 呈 指数 次 守 增 长 趋势 ， 运 
行 时 间 也 随 之 增加 。 阔 值 个 数 K=2，3 时 ， 各 算法 运行 时 间 
差距 较 小 。 当 取 到 较 高 阔 值 个 数 K=5 时 ， 各 算法 运行 时 间 差 
距 较 大 ， 其 中 ALO 算法 时 间 最 慢 。 当 K=3,4 时 PSO 算法 与 
ISCA 算法 运行 时 间 相 近 。 但 当 阔 值 个 数 提升 到 5 时， 本文 
算法 的 优势 得 以 体现 ， 其 运行 时 间 随 阔 值 个 数 增加 而 变化 的 


有] 


幅度 较 小 。 本 文 提出 的 ISCA 算法 ， 尺 管 随 着 闵 值 个 数 的 增 
加 运行 时 间 会 有 所 增加 , 但 在 病 值 个 数 K=5 时 仍 运算 时 间 优 (d)pig 
于 其 他 算法 ， 说 明 本 文 算法 效率 高 ， 能 够 较 快 地 完成 图 像 分 图 2 ”伯克利 原始 图 像 
割 任务 。 Fig.2 Original Berkeley images 
表 2 各 算法 最 佳 分 割 闽 值 
Table 2 Optimal threshold values obtained by each algorithm 
ISCA SAPSO SSA PSO ALO HSO 
图 像 K 
R G B R G B R G B R G B R G B R G B 
75 84 69 78 83 66 77 83 66 77 83 66 77 83 66 77 83 66 
139 50 1l9 146 149 125 145 149 125 145 149 125 45 49 125 46 48 125 
67 76 48 59 69 56 59 69 56 59 69 56 59 69 56 59 69 56 
3 104 11 70 05 05 85 105 105 85 105 105 85 05 05 85 05 06 85 
168 66 141 70 65 40 170 165 140 170 165 140 70 65 140 70 68 140 
53 58 49 53 60 51 53 60 51 53 60 51 53 60 51 51 61 51 
入 训 本 86 87 69 88 87 74 88 87 74 88 87 74 88 87 74 81 89 77 


80 60 142 174 158 136 174 158 136 174 158 136 74 58 136 74 57 137 


21 04 65 114 04 74 114 104 235 114 104 3 114 04 了 和 114 05 bg 


85 67 121 86 32 40 186 1 140 186 【72 140 86 72 140 86 3 141 

66 55 16 60 53 30 60 54 30 60 33 30 60 54 30 60 33 32 

09 93 66 06 93 65 106 94 65 106 93 65 06 94 65 06 94 67 

horseman 4 

77 50 103 58 39 99 158 140 99 158 139 99 58 40 99 60 41 102 

212 206 145 213 99 49 213 199 149 213 199 149 213 99 149 2153 204 L132 

39 41 22 52 42 32 $52 42 28 设 42 28 52 42 28 51 43 29 

2 2 35 84 72 66 84 72 59 84 但 59 84 Es 59 82 76 60 

3 ils 02 70 21 05 00 Wa 105 89 121 105 89 21 05 89 118 07 89 
55 34 132 66 46 31 166 146 123 166 146 123 66 46 123 63 49 123 

200 200 160 217 203 76 这 外 203 170 217 203 170 2 203 四 215 207 169 

2 00 06 用 04 04 67 104 104 67 104 104 67 04 04 67 04 03 67 


87 88 135 190 191 137 190 19] 137 190 191 137 90 NH 137 90 90 138 


210 200 149 210 204 3 210 204 153 1 204 153 210 204 153 209 204 153 


Tl 23 94 113 09 75 113 109 75 113 109 76 也 3 09 75 116 09 74 


bridge 56 43 141 58 Sl 116 158 i151 116 158 151 i117 58 31 116 61 村 116 
208 95 25 216 210 59 216 210 159 216 210 159 216 210 159 219 211 159 

63 79 3 66 64 30 66 64 2 66 64 23 66 64 2 69 65 31 

33 Il3 3 00 99 63 100 98 62 100 99 75 99 99 62 04 98 65 

5 47 29 85 34 32 99 134 131 09 134 132 116 33 32 99 36 33 101 

86 65 128 73 66 33 173 165 133 73 166 142 72 66 133 74 72 136 

231 216 157 222 216 66 222 216 166 222 216 170 222 216 166 224 215 164 

2 69 58 40 70 61 37 70 61 本 70 61 37 70 61 a 70 61 37 


45 1 了 7 EE 149 122 88 149 122 88 149 122 88 49 22 88 49 22 89 


67 44 93 162 132 124 162 132 27 162 131 97 62 32 27 62 32 99 


i 4 
PS 116 85 66 123 14 75 13 104 75 124 104 75 23 04 75 23 06 72 


89 py 108 180 148 114 180 148 114 181 148 114 80 48 114 80 51 109 


5 95 111 69 Sr 84 62 98 84 62 98 85 62 98 84 62 98 85 64 
117 55 106 136 115 90 137 115 90 137 116 90 37 ls 91 36 114 101 
184 218 160 188 156 126 189 156 126 189 158 126 89 Le 127 91 155 135 
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表 3 各 算法 运行 时 间 
Table3 Running time obtained by each algorithm 


图 像 时 间 (s) 
ISCA SAPSO SSA PSO ALO HSO 
2 3.0457 4.4085 5.0922 4.4982 7.6689 4.0567 
. 3 3.1123 4.7359 S77 4.5569 7.9659 5.6664 
a 4 3.5673 4.7697 5.2516 4.6020 8.8888 6.4161 
3 3.9151 5.3457 5.3464 4.7452 10.1538 7.2517 
2 3.3988 3.5215 4.0104 3.4799 6.4452 4.1924 
ee 3 3.6110 4.6704 5.2647 4.2391 7.6996 4.6256 
4 4.0773 4.9945 5.3161 4.6492 9.1019 4.8122 
5 4.1974 5.3839 5.8027 4.8519 10.8712 5.2971 
2 3.4149 4.0749 4.0320 4.6851 7.3420 4.3063 
biidge 3 3.6381 4.6919 5.4292 4.6422 8.4288 4.9381 
4 3.7298 4.1349 5.4154 4.3438 9.3038 5.1670 
5 4.2519 4.7859 5.4232 4.8032 10.6338 5.4593 
2 3.2110 3.8976 5.1616 4.4941 7.0998 4.3744 
. 3 3.6821 4.9248 5.2334 4.6945 7.9815 4.5586 
Ee 4 3.9859 4.9545 5.3763 4.7266 9.2002 4.7509 
3 4.1130 5.2493 5.7323 4.8389 9.8534 5.4511 
表 4 各 算法 分 割 图 像 指标 
Table 4 Indices of segmented image obtained by each algorithm 
PSNR SSIM 
图 像 K 
_ ISCA SAPSO SSA PSO ALO HSO ISCA SAPSO SSA PSO ALO HSO 
【en 2 4.3932 14.3737 4.3832 14.3832 4.3832 14.1392 0.3183 0.3169 0.3181 0.3181 0.3181 0.3087 
《 . 3 7.0466 17.0313 7.0313 17.0313 7.0313 16.7344 0.4638 0.4638 0.4638 0.4530 0.4638 0.4636 
wer 4 9.3209 19.1908 9.1908 19.1947 9.1908 19.2525 0.5737 0.5679 0.5679 0.5675 0.5679 0.5700 
有 21.6244 21.0907 20.9584 20.9790 21.0748 20.8383 0.6896 0.6522 0.6485 0.6494 0.6525 0.6398 
》 2 4.4799 14.3832 4.3832 14.3832 4.3832 14.3776 0.3252 0.3181 0.3181 0.3181 0.3181 0.3180 
3 7.0459 17.0313 7.0313 17.0313 7.0313 16.8233 0.4638 0.4557 0.4638 0.4638 0.4638 0.4498 
ea 4 9.7594 19.1917 9.1908 19.1878 9.1908 18.9160 0.5693 0.5675 0.5679 0.5679 0.5679 0.5498 
5 22.4564 20.9709 21.0748 21.0688 20.9601 20.7858 0.6874 0.6481 0.6525 0.6525 0.6485 0.6444 
2 4.2048 14.1361 14.1361 14.1361 14.1361 14.0884 0.4934 0.4918 0.4918 0.4918 0.4918 0.4889 
S . 3 6.7809 16.7379 6.7379 16.7379 6.7379 16.1889 0.6204 0.6191 0.6191 0.6191 0.6191 0.5887 
1 Pngee 4 9.3742 18.9590 8.9590 18.9588 8.9590 18.8898 0.7114 0.7113 0.7113 0.7070 0.7113 0.7082 
人- 3 20.3207 20.2488 20.2551 19.1329 20.2535 19.9317 0.7693 0.7680 0.7452 0.7637 0.7691 0.7593 
2 5.3413 15.3083 5.3083 15.3083 5.3083 15.3390 0.4216 0.4214 0.4214 0.4214 0.4214 0.4186 
. 3 8.1568 17.2892 7.2620 17.2335 7.2620 17.3977 0.5546 0.5246 0.5509 0.5510 0.5509 0.5346 
PE 4 9.5513 19.5423 9.5423 19.5413 9.5423 19.3491 0.6394 0.6377 0.6377 0.6375 0.6377 0.6354 
5 23.8224 21.3863 21.4014 21.4419 21.5022 21.9999 0.7016 0.6994 0.6993 0.6978 0.6996 0.6916 
多 闵 值 图 像 分 割 的 闵 值 个 数 越 多 ， 其 分 割 图 片 的 精度 越 。” 像 相 似 度 最 高 ， 说 明 该 算法 能 够 胜任 复杂 图 像 的 处 理 任务 ， 


高 ， 即 与 原 图 更 加 接近 。 为 定量 地 评估 各 优化 算法 在 分 割 效 。” 更 好 地 为 后 续 植 物 冠 层 图 像 的 识别 与 分 割 提供 支持 。 为 进 
果 方 面 的 能 力 ， 本 文 对 各 算法 的 PSNR 和 SSIM 指标 进行 比 ” 步 证 明 多 阔 值 分 割 方法 在 图 像 处 理 领域 中 的 优异 性 能 ， 能 很 
较 , 其 结果 如 表 4 所 示 。 首 先 分 析 峰 值 信 噪 比 PSNR, 当 K=2， 好 地 将 图 像 分 为 若干 个 具有 特殊 性 质 的 独立 区 域 ， 在 下 一 节 
3 时 ， 各 优化 算法 取 值 相近 ， 说 明 其 分 割 图 像 差 别 也 不 是 很 中 将 会 对 植物 冠 层 图 像 进行 实验 与 分 析 。 


大 , 但 相 比 之 下 本 文 算法 仍 有 一 定 优势 。 当 K=4, 5 时 , ISCA ”4.3 植物 冠 层 图 像 分 割 实验 

算法 的 分 割 指标 增加 明显 ， 对 于 tiger 和 horseman 这 种 目标 植物 的 冠 层 作 为 植物 最 先 接触 到 外 界 气 体 环 境 和 光照 的 
与 背景 相差 较 小 、 不 易 处 理 的 图 像 ， 数 值 有 明显 提升 。 说 明 ”部 位 ， 其 结构 形态 与 植物 的 生长 状况 密切 相关 。 为 提高 植物 
本 文 算法 不 仅 适用 于 多 维 函 数 极 值 求解 问题 ， 在 图 像 分 割 领 ” 生长 状态 评估 的 可 靠 性 ， 降 低 由 于 图 像 分 割 技术 不 成 熟 导致 
域 也 有 很 强 的 工程 实用 性 


E。 其 次 分 析 结 构 相 似 性 SSIM， 从 ”测算 误差 的 可 能 ， 本 节 对 植物 冠 层 图 像 的 分 割 进行 了 研究 。 
表 4 中 可 以 看 出 ， 随 着 阔 值 个 数 的 增加 ， 分 割 图 像 的 失真 程 ”选取 2 幅 植物 冠 层 图 像 进行 多 阔 值 分 割 实 验 ， 如 图 3 所 示 ， 
度 越 低 ， 更 加 接近 于 原 图 像 。 本 文 算法 得 到 的 指标 均 优 于 对 ”并 与 改进 分 水 岭 算法 以 及 模糊 聚 类 算法 进行 分 割 效 果 对 比 
比 算法 ，ALO 算法 随 着 阔 值 个 数 的 继续 增加 ， 其 结构 相似 性 ”P2-29。 为 得 到 较 好 的 分 割 效果 冰 值 个 数 以 及 聚 类 数 均 设 置 为 
指标 则 低 于 本 文 算 法 获得 指标 ， 说 明 ALO 算法 所 得 阔 值 与 ” 5。 实验 中 ISCA 算法 参数 的 设置 与 上 节 一 致 ， 模 糊 聚 类 算法 
理想 分 割 阀 值 相 比 仍 有 着 一 定 差距 ， 进 而 验证 了 ISCA 算法 ”和 迭代 终止 条 件 设 为 相 邻 两 次 迭代 差 值 小 于 0.001 。 
寻 优 精度 高 ， 能 够 得 到 较 好 的 分 割 闪 值 。 本 文 提 出 的 ISCA 为 从 分 割 图 像 中 分 割 出 植物 冠 层 部 分 ， 本 文 进行 了 如 下 
算法 与 其 他 对 比 算法 相 比 性 能 最 佳 ， 得 到 的 分 割 图 像 与 原 图 ”操作 : 


EM 
| 


到 


和 pa [一 了 叮 届 
( AX I\/ 二 4 上 口 十 | 
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Xl 


a) 从 分 割 得 到 的 K+1 个 区 域 中 选取 含有 植物 主体 的 
域 并 将 其 置 1， 其 余部 分 置 0。 

b) 对 选取 的 区 域 进 行 形态 学 处 理 ( 腐 蚀 , 膨胀 等 操作 )， 
在 得 到 的 图 像 中 提取 最 大 连通 区 。 上 述 操作 反复 进行 ， 直 到 
获得 外 轮廓 线 清晰 且 闭 合 的 图 像 。 

c) 对 得 到 的 最 大 连通 区 图 像 进 行 填 洞 处 理 并 与 原 图 像 矩 
阵 点 乘 ， 即 可 得 到 分 割 后 的 植物 图 像 如 图 4 所 示 。 

为 了 更 好 的 验证 本 文 提 出 的 算法 ， 将 最 大 业 作 为 评测 指 
标 对 比 各 算法 的 分 割 图 像 中 的 信息 量 ， 其 定义 如 下 : 
Entropy =—Ph logh—PlogP (13) 
其 中 : B 与 了 分别 为 图 像 前 景 、 背 景 像素 占 整 幅 图 像 的 概率 。 
炉 值 越 大 代表 图 像 所 包含 信息 越 丰富 ， 即 图 像 质量 越 高 。 
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(a)plantl (b)plant2 
图 3 植物 冠 层 图 像 
Fig.3 Plant canopy images 

各 算法 的 植物 冠 层 图 像 的 分 割 图 如 图 4 所 示 ， 从 图 4 本 
以 明显 地 看 出 ， 本 文 提出 的 基于 ISCA 算法 的 多 闵 值 图 像 分 
| 方法 优 于 其 他 分 割 算 法 ， 完 整地 将 植物 的 叶片 区 域 进行 分 
制 。 从 Plantl 分 割 图 像 中 可 以 看 出 ， 改 进 分 水 岭 算 法 分 割 效 
果 较 差 ， 存 在 欠 分 割 现象 ， 只 能 获得 部 分 叶片 。 模 糊 聚 类 算 
法 尽管 获得 叶子 大 致 轮廓 但 仍 存在 叶片 残缺 及 花 盆 区 域 ， 存 
在 过 分 割 现 象 。 而 相 比 之 下 本 文 算法 分 割 质 量 较 高 ， 能 够 将 
花 盆 与 冠 层 主体 很 好 进行 分 割 。 从 Plant2 植物 冠 层 的 分 割 结 
果 图 可 知 ， 分 水 岭 算 法 的 分 割 能 力 较 弱 ， 对 叶片 的 区 域 分 割 
不 完整 。 模 糊 聚 类 算法 能 够 分 割 冠 层 主体 区 域 ， 但 由 于 叶片 
的 受 光 不 均 导 致 分 割 效果 较 差 , 使 叶片 的 分 割 结果 存在 残缺 。 
本 文 算法 有 效 地 解决 了 受 光 不 均 题 ， 很 好 地 将 冠 层 主体 
归 为 一 类 ， 有 效 的 对 叶片 区 域 进行 了 分 割 |。 
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(a) 本 文 算法 (b) 改 进 分 水 岭 算 法 (c) 模 糊 聚 类 算法 


(d) 本 文 算法 (9) 改进 分 水 岭 算法 (DD) 模糊 聚 类 算法 
图 4 和 冠 层 分 割 图 像 
Fig.4 Segmented images of plant canopy 

各 算法 分 割 后 图 像 的 最 大 焙 值 和 运行 时 间 如 表 5 所 示 。 
从 表 5 中 的 最 大 炉 值 可 知 ， 改 进 的 分 水 岭 数 值 最 差 ， 本 文 提 
出 的 算法 优 其 他 对 比 算 法 ， 反 映 图 像 中 的 信息 更 多 ， 分 割 效 
果 更 好 。 从 CPU 运行 时 间 可 知 , 模糊 聚 类 算法 尽管 在 分 割 效 
果 上 优 于 改进 分 水 岭 算 法 ， 但 其 运行 时 间 较 长 。 而 本 文 算法 
的 程序 运行 时 间 最 短 , 对 于 K=5 这 种 高 阔 值 情况 其 运行 时 间 
较 短 ， 说 明 ISCA 算法 获得 分 割 图 像 的 速度 更 快 ， 实 时 性 更 
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好 ， 
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能 够 胜任 复杂 


能 够 在 保证 分 割 质量 的 基础 上 减少 运行 时 间 。 由 以 上 分 
可 知 ， 本 文 提 出 的 基于 ISCA 的 图 像 分 割 方法 性 能 较 好 ， 
图 像 的 分 割 。 

表 5 各 算法 指标 对 比 


Table 5 Comparison ofthe indices obtained by each algorithm 


图 像 


最 大 烂 运行 时 间 (s) 
本 文 算法 改进 分 水 岭 模糊 聚 类 ”本文 算法 ”改进 分 水 岭 、 模 糊 聚 类 


村 
BR 


lantl 0.55482 0.21022 0.43235 2.382716 76.933121 116.095853 


lant2 0.77895 0.22682 0.67370 4.568249 89.217856 199.667319 


5 


方 
从 


结束 语 

本 文 提 出 了 一 种 基于 改进 正 余弦 算法 的 多 阔 值 图 像 分 割 
法 ， 并 利用 该 算法 对 图 像 分 割 的 阔 值 选取 过 程 加 以 优化 。 
伯克利 图 像 的 实验 结果 可 知 ， 该 算法 有 效 地 提高 了 运行 效 


深 
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以 及 分 割 精度 ， 表 现 出 良好 地 适应 性 ， 很 好 地 平衡 了 探索 
开发 能 力 ， 为 图 像 的 进一步 处 理 黄 定 了 基础 。 从 植物 冠 


es 


守 叫 


] 
像 的 分 割 实验 可 知 ， 本 文 提 出 的 算法 能 够 完成 复杂 图 像 


分 割 。 今 后 将 会 寻找 更 多 有 效 的 智能 优化 算法 ， 并 尝试 与 混 
记 搜 索 策 略 结合 ， 将 其 应 用 到 多 闵 值 图 像 分 割 领域 以 提高 
的 寻 优 精度 。 
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